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Abstract. In the present paper we study the en- 
tanglement properties of thermal (a.k.a. Gibbs) states of 
quantum harmonic oscillator systems as functions of the 
Hamiltonian and the temperature. We prove the physical 
intuition that at sufficiently high temperatures the ther- 
mal state becomes fully separable and we deduce bounds 
on the critical temperature at which this happens. We 
show that the bound becomes tight for a wide class of 
Hamiltonians with sufficient translation symmetry. We 
find, that at the crossover the thermal energy is of the 
order of the energy of the strongest normal mode of the 
system and quantify the degree of entanglement below 
the critical temperature. Finally, we discuss the example 
of a ring topology in detail and compare our results with 
previous work in an entanglement-phase diagram. 



I. INTRODUCTION 

Entanglement [l], Q is one of the most characteristic fea- 
tures of quantum systems, and frequently the emergence 
of classicality is associated with the disappearance of en- 
tanglement from the state, i.e. separable states Q are 
deemed classical. This, and applications in quantum in- 
formation theory, e.g. [1, [H, [|| , have motivated a large 
literature on the problem how to distinguish separable 
from entangled states, see e.g. 0- 

Here we arc interested in the entanglement proper- 
ties of thermal (Gibbs) states of certain simple quan- 
tum systems: the starting point is the physically intu- 
itive observation that the ground state of a Hamiltonian 
of many interacting "particles" (we will adopt a more 
neutral, quantum information, language by calling them 
just "sites" ) is usually entangled, and so are the thermal 
states at sufficiently low finite temperature [|[ . For most 
well-behaved systems one would also expect the comple- 
mentary behaviour, namely that the thermal state at suf- 
ficiently high temperature loses all quantum character, 
and becomes separable, a concept formalised in [H . This 
is indeed true for finite systems of the spin lattice type, 
because their Hilbert space is finite dimensional, and in 
the limit of infinite temperature the thermal state con- 
verges to the maximally mixed state; then, results on the 



existence of a fully separable ball around the maximally 
mixed state can be invoked Q , to conclude that at large 
temperature the thermal state is separable. 

For infinite-dimensional or infinitely extended systems 
this type of argument is no longer available, yet the phys- 
ical intuition remains the same. In particular, in the 
present paper we focus on systems with a finite num- 
ber of coupled, infinite-dimensional harmonic oscillators. 
They provide a playground for testing ideas in quan- 
tum theory and quantum information. In particular, the 
class of quadratic Hamiltonians in all canonical coordi- 
nates {x\,p\, x 2 ,P2, ■ ■ ■) = (Ri, R2, ■ ■ ■ j Rin) provides a 
rich enough range of physical systems, yet the mathe- 
matics is simple enough so that almost everything can 
be understood exactly. 

Our objective is the thermal state of such a system and 
the question how entangled it is for a given temperature, 
especially the determination of the critical temperature 
at which all entanglement vanishes from the state. Ul- 
timately, we'd like to understand the continuum limit of 
this toy model: entanglement of quantum fields in the 
thermal state pH EH, EH , but the present simpler set- 
ting is interesting in itself, and will yield some decisive 
insights. 

Consider a system of n quantum harmonic oscillators, 
interacting with a quadratic self-adjoint Hamiltonian, the 
most general form of which is 



2n 



H 



X] C 3k R j R k + <l.l>',- 
j,k=l 



2n 



(i) 



with coefficients Cjk — c* k - for the quadratic terms and 
displacement coefficients dj = d* for the linear terms. 
This Hamiltonian might, for example, describe a system 
of interacting oscillators in a lattice with some internal 
interactions; for the moment we will however look at the 
most general case. 

As for spin systems, the thermal states at low temper- 
ature T, 



-/3H 



with Z = tr[e 



(2) 
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arc in general entangled [H], where (3 = (fc^T) -1 is 
the inverse temperature. Conversely, it is intuitively ex- 
pected that for sufficiently high temperature the thermal 
states are separable. In this paper we prove this intu- 
ition mathematically, and give quantitative bounds on 
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the critical temperature at which separability kicks in, in 
Section IIIII We proceed to discuss the special case of a 
translation-invariant Hamiltonian (e.g. a ring of coupled 
harmonic oscillators) in Section IIVI for which the criti- 
cal temperature bound is shown to be tight. We provide 
an in-depth discussion of these results when applied to 
the harmonic ring with nearest-neighbour couplings in 
Section [V] In Section IVII we give a brief discussion of a 
simple entanglement measure which can be evaluated ef- 
ficiently for thermal states at all temperatures and relate 
it to an established measure, the (Gaussian) entangle- 
ment of formation. In the final Section IVlII we close with 
a discussion of the physical meaning of our results and a 
comparison to previous work. 



II. CONTINUOUS VARIABLE STATES 

Quadratic Hamiltonians are the working horse of theo- 
retical physics and many a theory is based on the under- 
standing of these simplest of models. So it is in the case 
of continuous variable states, which have been studied ex- 
tensively on the model of Gaussian states (for reference, 
see for instance [2 EB|)- Gaussian states can in par- 
ticular be represented as the thermal states of a system 
with a quadratic Hamiltonian |T]) and we can therefore 
use the well-developed machinery of analysing Gaussian 
continuous variable states for our purpose. In partic- 
ular, the thermal states under discussion are uniquely 
specified by the expectation values of all position and 
momentum variables, i.e. (Rj) which vanish for all j, 
and by the quadratures (covariances) (h{Rj, Rk}) = 
(i(i?ji?fc + RkRj)) = 7jfe. The latter form a 2n x 2n- 
matrix, 7, which is real symmetric and satisfies the so- 
called uncertainty relations 7 + ir-cr > 0, with the sym- 
plectic form a determined through the commutation re- 
lations [Rj,Rk] = ihajk and given explicitly by 
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1 

-1 



(3) 



These uncertainty relations are sufficient for a covariancc 
matrix (CM) to define a bona-fide quantum state and 
the first moments can be changed arbitrarily by applying 
local, i.e. single-oscillator, Weyl-displacement operators. 
Therefore, only the quadratic part of the Hamiltonian 
in Eq. ((T|) is relevant and the coefficients dj can be set 
to zero. Thus the entanglement properties of pp are in 
general determined by the covariance matrix 7 and also 
higher moments; in the Gaussian case by the covariance 
matrix alone. 

Remark 1 Even though the following will be obvious 
to some readers, let us here make a remark regarding 
the choice of the basis for the class of systems under 
discussion to avoid a common confusion. The original 
physical oscillators described by the canonical observ- 



sight in an alternative mathematical description of the 
system. In fact, any symplectic transformation S trans- 
forms into a collection of n abstract modes, R lT = SR T , 
that will be coupled in a different way from the original 
oscillators. This change of coupling results in a change 
of the entanglement properties between the underlying 
modes; for instance the normal modes are by definition 
completely decoupled and therefore never entangled at 
any temperature, see Fig. [TJ It is hence useless to speak 
about complete separability in an abstract way. There 
is no such thing - separability may be had in one basis 
or decomposition but not in a rotated, different basis at 
the same temperature, depending on the choice of the 
division. It is not only a grouping issue of how many 
oscillators one puts in several groups - the issue occurs 
already when the choice of the representation is made - 
real space (physical oscillators) or k-space (normal modes 
or phononic modes, compare [l6|) or any "intermediate" 
set of modes. 



III. HIGH TEMPERATURE IMPLIES 
SEPARABILITY 

We will show that for sufficiently large T, i.e. sufficiently 
small (3, the thermal state is separable in the sense that 
it is a mixture of single-oscillator Gaussian states. This 
holds true for any arbitrary choice of basis of the system, 
i.e. separable with respect to the original spatial basis 
of physical oscillators (i?) or w.r.t. the normal basis of 
decoupled oscillators (R 1 ) or any other basis that can be 
reached by symplectic transformations. 

Separability of a Gaussian state can be formulated in 
terms of the covariance matrix alone. Werner and Wolf 
[T^ have shown that a bipartite Gaussian state, with 
its covariance matrix denoted 7, is separable if and only 
if there exist Gaussian states on the local systems with 
covariance matrices r\A and t\b such that 



7 > VA © Vb- 



(4) 



Note here, that matrix inequalities for matrices con- 
taining physical units, such as position and momentum, 
should be read as being evaluated with respect to vectors 
\v) which themselves carry units; for instance 7 > r)A®rjB 
shall be read as (i>|7|t>) > (v \t]a © Vb\v) for all vectors 
\v). 

In easy generalisation of Eq. ([5]), our thermal state pp 
of n oscillators with covariance matrix 7^ is fully sepa- 
rable into single physical oscillators if and only if there 
exist Gaussian states of the individual oscillators, with 
covariance matrices r]j (j = 1, . . . , n) such that 



7/3 



> 



(5) 



ablcs Xj 



i?2j-i and pj ~ i?2j may disappear from 



We shall show that for any arbitrary n-tuple of one-mode 
covariance matrices rjj (each a 2 x 2 matrix) this even- 
tually becomes true for high enough temperatures. The 
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fully separable into thermal normal modes 
fully separable into some "intermediate" modes 

fully separable into physical oscillators in space 



FIG. 1: Separability is a concept depending on the choice of 
the basis of the system. Shown here is the set of all states 
(black set) which a system of n harmonic oscillators can re- 
alise in principle. We wish to look at only the thermal states 
for a given Hamiltonian and vary the temperature T (black 
line). On a journey from high temperatures towards low tem- 
peratures we pass through the following stages: At extremely 
high temperatures we expect that the system is fully separable 
for all possible decompositions into n modes. This is because 
high temperature provides so much energy that the oscillators 
lose all correlations and behave completely independent from 
each other in any representation. When the temperature is 
lowered the system reaches a point where something in the 
chain must become entangled. The real space representation 
referring to the actual lattice sites turns out to be the repre- 
sentation where entanglement can be found first, at least for 
Hamiltonians with translational invariance of the coupling (as 
discussed in Section ITV|| which implies that the derived bound 
for the critical temperature T cr it, see Eq. (|43[1 in Section Hvl is 
in fact tight. In other words, for temperatures below T cr it the 
thermal state picks up entanglement between sites in space 
and leaves the set of fully separable sites (red). However, any 
"intermediate" mode representation between physical-space 
and the normal mode decomposition will pick up entangle- 
ment only at smaller and smaller temperatures (blue). The 
extreme case is of course the normal mode basis for which no 
temperature is low enough to produce entanglement and any 
thermal state remains separable (green). 



reasoning is essentially as follows: The matrix ©™ =1 Vj 
has a fixed largest eigenvalue, denoted A max . However the 
smallest eigenvalue of turns out to diverge as (3 — ► 
and therefore the matrix inequality must hold above some 
finite critical temperature T cr ; t . 



Theorem 1 For every quadratic Hamiltonian of the 
form J3p, there exists a temperature T cr j t such that for 
all T > T cr it the thermal state pp with (3 = (^T) -1 is 
fully separable. 



Proof: There exists a symplectic transformation ma- 
trix S of the canonical variables R' T = SR T , which diag- 
onalises the Hamiltonian (JTJ) with respect to the primed 
variables, 



n n 1 



a 2 j Rfj 



(6) 
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where S is chosen such that J2 m n ^km 
Any thermal state pp can thus be expressed as a (tensor) 
product of Gaussian states of single harmonic modes in 
the corresponding normal basis, 



-fSHj 



(7) 



The condition that the new variables are canonical is 
expressed by the symplectic condition S a S T = a and 
the covariance matrix in these new coordinates is simply 
7^ = S7/3 S T . Since S is invertible, the separability con- 
dition can be brought over to the primed matrices. 
The reason is that even though the symplectic transfor- 
mation S may change the eigenvalues of a positive matrix 
it will not change its positivity. Therefore the separabil- 
ity condition ((S]) becomes 



S~/pS T > s 




To, 



(8) 



for some constant covariance matrix Tq with largest 
eigenvalue Ao- On the left hand side, however, we have 
the covariance matrix of a tensor product state, i.e. the 
covariance matrix 7^ is the direct sum of 2 x 2-blocks, 



(9) 



where 7j(/3) are the covariance matrices of the thermal 
states of the harmonic modes each with their Hamiltoni- 
ans Hj, at inverse temperature [3. These blocks are given 
explicitly by 



■ coth 



(3hwj 





rrij u)j 

rrijUJj 



(10) 



where u)j = ^/o^j-i ot2j is the frequency of each of the 
modes and rru = 1/ a.2j is the mass of the physical oscil- 
lators (see [18| for reference). 



Introducing the length unit Lj we rescale the variance 
in position Aj — 2 m u- » ana * w ith h/Lj being the ap- 
propriate fundamental momentum unit the variance in 



momentum becomes £>,■ = 



hi 



The final unit-free 



3 2 

covariance matrices for each independent normal mode 
are then 



Aj coth 



(3hujj 





Bj coth 



(11) 
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with the unit-free eigenvalues on the diagonal which in- 
deed go to infinity for [3 — > for any arbitrary finite 
unit-length Lj. Alternatively, one can calculate the 
symplectic eigenvalues of the unit-free matrices as 



coth 



2 ' 



so that s., > i for all 



Sj = yjdet 7j- (/?) - 

temperatures and Sj — > oo for j3 — > as well. In other 
words, for suffiently small (3 (i.e., sufficiently large T), 
the matrix inequality (|5|) will be satisfied, proving the 
thermal state separable. □ . 

So far we could confirm the intuition that for (3 — > 
the thermal states p/3 take on a separable configuration 
for any Hamiltonian and irrespective of the chosen mode- 
decomposition. However, we want to improve our argu- 
ment such that we can predict the quantitative critical 
temperature for a particular (class of) Hamiltonian at 
which the thermal state undergoes the transition from en- 
tangled to separable. For this purpose, let us go through 
the above proof once more, but now paying attention to 
the freedoms we have there to get as small as possible 
a bound on the critical temperature - thus turning the 
quest for the exact value into an optimisation problem. 

Let S be again the symplectic transformation that di- 
agonalises 7/3, and rescale to the unit-free version of the 
matrix, 



S 1 units ~ 
7/3^7/3' 











(12) 



Now consider all unit-free covariance matrices Tq that 
can be reached by first choosing individual covariance 
matrices for all modes [the m for j = 1 , . . . , n in 
Eq. (|SJ)], then, secondly, transforming them with S, 

Tq = S f©" =1 r/jj S T , and finally rescaling to the unit- 
free matrix Tq. We require that every eigenvalue of 
j'g [see Eq. (|12j) ] must be larger than every eigenvalue 

of To, in particular larger than the maximal eigenvalue 
Ao = A 

max (To), 



1 , Ptwi 
s.j = — coth > Ao 



for all j = 1, 



(13) 



Therefore the separability condition for the thermal 
states pp is sufficiently satisfied when (3 obeys 



&rit < 



1 



fojJr, 



■In 



2A + 1 
2A -1 



(14) 



This still is a very rough bound for the critical temper- 



IV. SYMMETRIC SYSTEMS 

The recipe of the previous section turns out to be 
pretty good for systems with a lot of symmetry, and some 
additional structure of the Hamiltonian: Following Audc- 
naert et al. (20| , we consider systems with a Hamiltonian 
of the form (now in position and momentum coordinates) 



H 



^2m^ 

3 



Em 

jk 



(15) 



with Vjk = Vkj = Kr(j) „fjA f° r an Y permutation ir. (A 
particular such system is the translation invariant ring 
of n equal harmonic oscillators, coupled via a nearest- 
neighbour interaction u> and exposed to on-site potential 
traps with strength S, see Section W\) 

Note, that the Xj all commute, hence we may assume 
that the potential matrix is real symmetric and for 
physical reasons it needs to be positive semi-definite so 
that the energy is bounded from below. Then it is shown 
in [20( | that the covariance matrix of the thermal state 
at inverse temperature (3 can be written as a direct sum 
of a position part, 7x(/3), and a momentum part, 7p(/3), 
7/3 = 7x(#) ®1p(/3), with 



ft „ 1/2 , pnv 1 / 2 
7*(/3) = —F- 1/2 coth-" 



2m 
<nh 

T 



7p(/?) = ^ 1/2 coth 



2 

(3KV 1 / 2 



(16) 



Writing the covariance matrix in this way requires re- 
ordering of the canonical variables into all the positions 
first followed by all their corresponding momenta. The 
observation in [20 ] was, that because of the block nature 
of 7/3 it can actually be diagonalised by an orthogonal 
symplectic transformation of the special form S®S, with 

xp 

an orthogonal n x n-matrix S. As a consequence of sym- 
plecticity, S ffi S decouples the oscillators in the sense 

xp 

that in the new canonical coordinates x'j = SjkXk 
and p'j = SjkPk, the Hamiltonian assumes the form 



H 



E 



2m 



E 

3 



2 U 3 X 3 



(17) 



In the canonical coordinates . 



p'j the covariance matrix 



70 becomes a direct sum of blocks (compare (jTUJ) ) 



2rnuj- 



coth 



2 



coth 



(3hujj 
2 



(18) 



ature, T cr it = 



as there is no reason in general to whcre the masgcs arc ^ thc same and the frequencies 

jj arc obtained by diagonalising y/V ; indeed, the us 2 are 



expect that we have to compare thc smallest eigenvalue 
on the left with the largest on the right, i.e. this condi- 
tion may not be necessary for separability. Nevertheless, 
a slight refinement of this approach can give reasonable 
and intuitively appealing bounds on the critical tempera- 
ture, which arc in fact sometimes exact, as we shall show 
in thc next section. 



just the eigenvalues of the real symmetric matrix V . 

We intend to compare 7/3 with (B" =1 Vo, a direct sum of 
2 x 2-blocks 770, which we argue can be chosen w.l.o.g. as 



'/() 







2mujQ 

q mhujp 



(19) 
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This is straightforward since first, the matrix ordering is 
preserved under the map that sets all "mixed term" co- 
variances to zero. Secondly, the state can be assumed as a 
minimal uncertainty state - otherwise we could decrease 
either of the diagonal entries of t/q. As a consequence 
of orthogonality of S © S, 0" = x is of the exact same 



;f ■ p 

form in the transformed coordinates x'j, p'y 



S (B S 

xp 



S (B S 

xp 



= ®4» 

.7=1 



(20) 



with tj'q having a possibly re-scaled frequency co' , but 
since lvq is still arbitrary we drop the prime immediately. 

This implies that we get 2n inequalities which have to 
be satisfied by /3 to guarantee separability: for all j, 



coth 



> — and 

UJQ 



coth 



> 



u 



-■ (21) 



(Note, that by coth a; > 1 only one of the two inequal- 
ities expresses a nontrivial constraint for every j, so we 
have effectively only n inequalities.) A brief calculation 
(distinguishing the cases of ujj < u>q and ujj > u>o) shows 
these constraints require 



jj < min 

3 



1 



in 



UJj + LUQ 



UJj - LJ Q 



Introducing the scaling function 



s(x) 



we can rephrase this as 



1 



■In 



1 



1 . / Wi 

p < mm s — 



(22) 



(23) 



(24) 



The good thing here is that we see how the critical tem- 
perature depends only on a few universal parameters and 
items: essentially, j3 is given by j^-, with the scaling 
function giving a constant multiple. The properties of 
s are very simple: it is monotonically increasing in the 
interval [0; 1), and decreasing in the interval (l;oo), and 
s(x) > 2 for < x < 1 (see Fig. 0). Indeed, the Taylor 
expansion of s at x = is 



s(x) = 2J2 



— x 2k 



k=0 



2k + 1 



and it obeys the functional relation 



s{x) 



1 



(25) 



(26) 



This implies that actually only the extremal values out 
of the cigcnfrequcncy spectrum of H . {ujj}, become rele- 
vant: 



a , s 1 
p(u>o) = - — mm 



L0 



(27) 
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FIG. 2: The scaling function s(x) is monotonically increasing 
left of the singularity at x = 1 and decreasing right of it. 
This function encodes the multiplicative offset from -^huio 
for the critical temperature. The dashed blue line indicates 
the cutting point x = 1.199678 ~ 1.2 where s(x) passes below 
its initial value of 2. 



where now /3(wo) is an upper bound so that f3 < (3{uio) im- 
plies separability. Because we want to use this to find (or 
bound) the minimum temperature at which the thermal 
state is separable, we would like to maximise (3{uiq) by 
choosing loq appropriately. For example, we find a physi- 
cally intuitive bound for the critical temperature at which 
all entanglement vanishes from the oscillator system: by 

choosing uj = w max , we find /3(w raax ) = 7— ^ s 

with values of s only in the interval [0; 1] where they are 
always larger or equal to 2. Hence, if the thermal energy 
available to the system, ksT, is at least half as big as the 
largest energy of an eigenmode of the system, /iw max , the 
system will be separable. 

By a simple variational argument it is seen that to 
maximise 0(loq) it is necessary that s (^73^ = s (^^^^j , 
and since this has a unique solution for lj , this value 
is indeed the one attaining the maximum; i.e., if j3 < 



max^g P(uiq), then the thermal state at inverse 



temperature (3 is fully separable. In fact, the solution 
to the equation for s depends only on the ratio of the 
spectrum, r := ^ max > 1, and we can summarise our 
findings so far in the following theorem. 



Theorem 2 For all inverse temperatures /3 < /3 cr it := 
j-j — <r(r), the thermal state is separable. Here, cr(r) = 
ts{t) with the unique 1 <t<r such that s(t) = s (-). 



Proof: We simply go with u> 



into the above 



argument; there, it is straightforward to see that nothing 
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is gained by choosing luq outside the interval [w m i n ; w max ] 
□ . 

By looking at the graph of s once more (Fig. [2]), we find 
that <j{r) is monotonically decreasing with r. Its infimum 
is attained at r = +00, where t > 1 is the unique solution 
to the equation s(t) = 2, yielding t = 1.199678 . . .; hence 
cr(+oo) = 2 ■ 1.199678. . . = 2.399357. . .. On the other 
hand, for very small r (i.e., close to 1), we can equally 
easily understand the asymptotics of o~(r): by choosing 
t = ^[r - a good choice as can be seen by looking at the 
functional equation (|26p - we find 



cr(r) > min 




for r = 1 + e. Indeed this is asymptotically exact as 
r -> 1 + 0; see Fig. H 



r) _ 
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"V^^^ lower bound 


1 : 


asymptotics for r — > 1 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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the thermal state at inverse temperature (3 > /3 cr j t is not 
fully separable into physical oscillators. 

Proof: For this, all we have to show is that the ansatz 
of comparing 7 to (Bj =1 rjo is already the most general 
we need to consider. Indeed, the covariance matrix 7 
inherits the symmetry of now on the level of pairs 
of canonical coordinates; but then, if 7 > Tjj, also 

7 > W^ti) f° r an y P ermu t a ti° n f S G, and since G 

acts transitively, we find by averaging that 7 > ©™ =1 ?70j 
with 

1 1 - 

Hence, all of our calculations above were indeed 
w.l.o.g. and there is no way we could have found a smaller 
critical temperature □ . 

Examples of transitively acting symmetry group of the 
Hamiltonian abound: for instance the shift-invariant har- 
monic ring considered in the following Section [Vj or any 
other periodic lattice in higher dimension. 




0.5 1 r 



FIG. 3: A plot of a as a function of r comparing it with its 
two asymptotic lower bounds: The constant 399357 . . . ~ 2.4 
in the limit r — * 00 (see also Fig. [2} and the expression from 
Eq. ([28j) for r -> 1. 

While all this might feel quite ad hoc (though giving 
a reasonable bound), we want to show now that the de- 
scribed method to obtain /? C rit actually yields the exact 
cutoff point for all systems of the type (fT5|) with suffi- 
cient translation symmetry. To be precise, let G be the 
group of site permutations which leave the Hamiltonian 
invariant. Then we have the following complement to 
Theorem H 

Theorem 3 Under the previous assumptions on the 
Hamiltonian, and if G acts transitively on the sites of 
our n oscillators ( meaning that for every two sites j and 
k there is a permutation tt G G such that 7r(j) = k), then 



FIG. 4: This is the universal entanglement phase diagram of 
translation symmetric systems of type Eq. (|15|l : The bound- 
ary marks the line between the system being fully separable 
into separate sites for high temperatures, and it being (some- 
how) entangled for low temperatures (green shaded). The 
temperature on the vertical axis is scaled with fiw milx /<;B, the 
horizontal axis is the inverse of the spectral ratio, i = ""*'" . 

This result shows that the critical temperature for the 
class of systems considered in the theorem depends really 
only on the top end of the normal mode spectrum w max , 
and the ratio r of the spectrum; in Fig. [5] we interpret 
this result as describing a universal phase diagram of the 
presence of entanglement in the system. 

Remark 2 Here we point out why all the assumptions 
of Theorem [3] are necessary in our proof: we have to 
impose translation symmetry to be able to reduce to the 
case of checking the separability criterion on a direct sum 
of identical blocks 770. Then we need the assumption 



7 



of type (fT"5|) Hamiltonian for two steps: the first is to 
be able to get rid of the off-diagonal entropies of rj , so 
that © J= i f]o is the direct sum of two multiples of the 
identity (collecting the position and momentum variables 
in contiguous blocks); the second is to have an orthogonal 
symplectic transformation of the form S(BS diagonalising 

xp 

H, which leaves ?yo alone. 



V. THE HARMONIC RING 

Let us take a particularly symmetric system, the trans- 
lation invariant ring of n equal harmonic oscillators, cou- 
pled via a nearest-neighbour interaction u> and exposed to 
on-site potential traps with strength 5. The Hamiltonian 
of such a quantum system is 



H 




(30) 



where m is the mass of each physical oscillator. (The 
convention is that indices are understood modn.) Diag- 
onalising via discrete Fourier-transformation leads to a 
Hamiltonian with n decoupled oscillators with the eigen- 
frequencies 



2u\ 



' sm 



2 



2w 



withj = 0, . . . ,n- 1. (31) 



These eigenfrcquencies are on the scale of the interaction 
frequency u> but depend heavily on the trapping potential 
5 as well. In the special case 5 = 0, a particular spatial 
symmetry is assumed, i.e. the coordinates of the centre of 
mass of all oscillators remains completely free. However, 
for any finite, ever so tiny, S this degree of freedom is 
harmonic. The scale of energy in this discussion is given 
by the gauge frequency loq where w m i n < < w max 
can be chosen in an optimal way to achieve the tightest 
bound. Roughly, ujq will be u) max , and the correction to 
that depends only on the ratio of the spectrum r = " max 
which itself depends on S. 



According to Eq. ([31]) the minimum frequency which 
gives the cost of the lowest oscillation is w m in = S and 

the maximum frequency is u) max = 2uo\ \ 



4- . The 



ratio between these two is then r = 
critical temperature becomes, by Theorem [3l 



1 + (2".y and the 



fc B T r 



B -tcrit 



cr(r) 



H6- 



a(r) 



which is an implicit function of — 
5 — > 0, the ratio is approximately r 
therefore 



via r = r 

5 



(32) 

§-). For 

3> 1, and 



l 

L2' 



(33) 



entanglement 
btw. two neigh- 
bouring blocks, 
of size n/2 




5 1 LP 

21nS/w 



ML 



completely 
separable into 
independent 
sites 



temperature 
thresholds 

for 5 = 



ML 



L2k B 



FIG. 5: Schematic of how the entanglement is distributed be- 
tween the spatial sites for varying S (n is thought to be large). 
When 8 = 0, bipartite entanglement resides between near- 
est neighbours (green) for temperatures up to T n . n . = jT^" 
and to a lesser degree and for smaller temperatures between 
non-nearest neighbours [l(|. On the other hand, according 
to Eq. (|33[) , above T cr i t = hiv/1.2kB there is no entangle- 
ment for any possible split (yellow). Between T n . n . and T cr it 
entanglement must be present but not simply between two 
single sites. Audenaert et al. [2(| have reported numerical 
results showing that entanglement above T n . n . exists between 
two neighbouring blocks of size n/2 up to a temperature of 
^blocks ~ 2.8fttj/\/20 fcs for the case u> = \/20<5 (the unit- 
less value 2.8 appears to be the critical value according to 
Figure 12 in [201 ] : our expression is obtained by retaining all 
units). When S is large compared to ui, the frequency spec- 
trum Eq. (|31|l shrinks to essentially one dominating frequency, 
ujj w 8. The mixing of frequencies occurring at increasing 
temperatures becomes therefore irrelevant and entanglement 
persist up to higher temperatures, the higher 8 becomes. By 
Eq. (|34[) . the asymptotic behaviour is B R a 



5/u 



2 1n(S/u) ' 



where <t(2j) r* 2 ■ 1.2 for small 5. In the limit where 

6 — > oo, the ratio becomes r ~ 1 + and a if) ~ In ^ 
and therefore 



k B T, 



B -L crit 



5/u 



21n(r5/w)' 



(34) 



It turns out that the order of magnitude of tempera- 
tures where entanglement can occur is indeed set already 
by the nearest-neighbour entanglement, with full separa- 
bility kicking in at only roughly double the temperature 
where nearest neighbour (n.n.) entanglement disappears. 
A qualitative picture of this behaviour in comparison to 
the nearest neighbour entanglement is shown in Fig. 



VI. MEASURING THE ENTANGLEMENT AT 
SUB-CRITICAL TEMPERATURE 

Now that we know a pretty good bound on the 
entanglement-critical temperature for general systems, 
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and indeed the exact value for translation-symmetric sys- 
tems of the type (fT5|) , we proceed to an attempt to quan- 
tify the amount of entanglement below the critical tem- 
perature. In this section, we stick to systems of the latter 
type, even though one obtains bounds on the entangle- 
ment for more general systems in the same way as before. 

The most straightforward thing to do is to measure by 
how much we fail to satisfy the inequality 7 > ?7j, 
compare Eq. ^j. E.g. one may set 

n 

P(7) := max{0 <p<l:7>p0i|j for CMs 

3=1 

(35) 

and consider — logP(V) as a measure how far 7 is from 
being separable (see [21(, where the above is considered 
for bipartite splits and denoted V , and [U[). In [2l[ it is 
shown that this is indeed a Gaussian entanglement mono- 
tone (i.e. a measure monotonic under Gaussian LOCC). 

As before, by translation symmetry and the absence of 
mixed terms on the covariance matrix, we may assume 
w.l.o.g. that all r)j in Eq. (|35[) equal some 



q m hujo 



(36) 



This means that, with r max denoting the maximum of 
the r k , E(Y) = J2k E (y( T k)) > £( 7 '(w)), while 

pA\B {l) > pA\ B{r/l) = min pA|i3 ( y (Tfc)) 

k (40) 
= P A|B (V(w)) =e— , 

where the first inequality follows from 7 > 7' and the 
definition of P, and the first and third equality sign are 
facts proven in [U. Hence, -\nP A \ B (j) < r max . 

Using the formula for the entanglement entropy of 
the two- mode squeezed state j'(t) [24], we arrive at the 
bound 



^ |B (7)>M-lnP A|B (7)), 



(41) 



with 7i(r) := cosh 2 r log cosh 2 t — sinh 2 rlogsinh 2 t, a 
kind of "hyperbolic entropy" . This results - after some 
elementary transformation - in a relation 

i?G |S (7)>-21ogP^ s (7) + A(P^ s (7)), (42) 

where A is a monotonic and convex correction (in P) with 
the properties A(l) =0, A(0) = ^-2 = -0.557304.... 



VII. DISCUSSION AND CONCLUSIONS 



yielding 



P(7) = max min 

u>o j — 1 ... n 



3 , — coth 



UJ-, 

— coth 



u 



f3hujj 



U) 



(37) 

This is not as easy to evaluate as the critical tempera- 
ture (where P(7) becomes 1) but still remarkably simple. 
Note that as before for the critical temperature, the value 
depends only on the frequency spectrum of H, not on the 
details of the phononic modes. Now, however, all the fre- 
quencies play a role in P(7) as a function of the inverse 
temperature /?, not only the smallest and the largest. 

As a side remark, if we do the optimisation of Eq. ([33)) 
for a bipartite cut A\B of the system, P A \ B is related to 
other nice entanglement measures, specifically the Gaus- 

A I B \ — ^1 

sian entanglement of formation E G (7) [22|, for covari- 
ance matrices 7. Namely, as shown there, 



E. 



A\B 



(7) = inf {E(j) : 7 > 7'and7'purc} , (38) 

7' 



where £(7') is the entropy of entanglement of the pure 
entangled Gaussian state with covariance matrix 7'. Fur- 
thermore it is known [22j that a pure Gaussian state 7' 
is equivalent (via Gaussian local operations) to a direct 
sum of two-mode squeezed states, 7' = ® fc "f'ijk), where 



7'(r) = 



coshr 


sinhr 






coshr 


— sinh r 



sinh t 

— sinh t 

cosh r 
cosh t 



(39) 



We have shown that in all systems of finitely many 
(harmonically) coupled harmonic oscillators, the thermal 
state at sufficiently high temperature is fully separable. 
We note, that our argument was quite independent of 
the actual Hamiltonian, and may still hold for a much 
wider range of systems. For systems whose Hamiltonian 
has separate kinetic and potential terms for the phys- 
ical oscillators in space we have found a bound on the 
corresponding critical temperature, 



^'sTcrit < 



a(r) 



(43) 



where a(r) is a universal scaling function of the spectral 
ratio r = ^ max between the maximal (minimal) eigenfre- 
quencies of the Hamiltonian w max (w m i n ). This is a quite 
intuitive result, as it says that the temperature has to be 
just high enough so that all the normal modes are partic- 
ipating sufficiently in the thermal state. In general one 
can of course not expect the above formula to give the 
correct cutoff temperature for entanglement, as it neither 
refers to the full eigenfrcquency spectrum, nor at all to 
the actual form of those eigenmodes themselves. For ex- 
ample, as we pointed out in Section|TTl a system of decou- 
pled oscillators with eigenfrequencics ojj is disentangled 
at all temperatures. For Hamiltonians with sufficient 
translation symmetry, however, such as periodic chains 
or lattices, the right hand side of Eq. (|4*3")l turns out to 
give exactly the critical temperature. The intuition why 
this can happen is that the translation symmetry of the 
Hamiltonian sufficiently determines the entangled struc- 
ture of the eigenmodes. 



9 



For the periodic harmonic chain, our results made 
for some interesting comparison with earlier work, such 
as and [TB]- There, only two parameters govern the 
behaviour of the system, the interaction strength ui and 
the on-site potential S, see Eq. (|30|) . We find a phase di- 
agram of different entanglement/separability regimes as 
shown in Fig. \5\ The general case of quadratic Hamil- 
tonian is qualitatively similar in terms of the maximal 
cigenfrequency w max and the ratio r = ^ max . 

Based on the Gaussian separability criterion, we were 
then able to even quantify the entanglement at sub- 
critical temperatures via a certain distance to the sep- 
arable set. The most amazing aspect of the resulting 
formulae, Eqs. (|57|) and is that these expressions 

are given in terms of the normal mode spectrum alone. 
This is because the entanglement properties of the nor- 
mal modes themselves are sufficiently captured by the 
translation symmetry It should be noted, however, that 
the P-measure of entanglement, though related to the 
(Gaussian) entanglement of formation for bipartite cuts, 
is not an extensive quantity, hence it can not be used to 
extend the results of [24[ on entanglement-area laws for 
the ground state of harmonic systems. 

We had to leave open a number of intriguing questions, 
among them the calculation or at least estimation of an 
extensive entanglement measure, such as (Gaussian) en- 
tanglement of formation between two disjoint subsets of 
oscillators, especially the case of a region and its comple- 
ment in a lattice, as a function of the temperature and 
the shape of the regions, for example like the area laws 
for ground states discussed in [24|]. One can also spec- 
ulate about the existence of bound entanglement in the 
thermal state, more precisely, temperatures at which the 
state is PPT (positive partial transpose) and entangled. 
Related to this is the more general question about the 
kind of entanglement in the thermal state, in particu- 
lar whether it can be characterised by bipartite entan- 
glement or if it is genuinely multipartite. A recent paper 
[25( addresses exactly this issue and our results combined 
with those in [25[ prove, that fully PPT entangled states, 
i.e. entangled states such that all bipartite partitions are 
PPT, cannot be obtained for the harmonic systems inves- 
tigated here. This is because the critical temperatures for 
the negativity of the even-odd partition in [251 ] coincide 
with our Eq. (f27|) for full separability. Thus, the corre- 
sponding thermal states become PPT and fully separable 
for the same temperature. 

Another, deeper and physically more interesting, ques- 
tion is whether the regions in our entanglement phase di- 
agram, Fig. have any actual interpretation as different 
phases of the system in the sense of physically radically 
different orders. Clearly, we do not expect all the differ- 
ent boundaries to be phase transitions crossovers, just as 
we don't expect all possible entanglement measures to be 
thcrmodynamically relevant . 

For certain approaches to the quantum-to-classical 
transition stipulating its being connected to the loss of 



entanglement the following observation may be relevant. 
For a macroscopically large system, we would expect that 
the quantum description might have faster and faster nor- 
mal modes, requiring that the temperature guaranteeing 
separability goes to infinity. This clearly means that for 
a reasonable emergence of classicality, full separability of 
the state is asking too much; on the other hand, it can be 
interpreted as saying that even large and hot systems re- 
tain some quantum character hidden in very high energy 
degrees of freedom. 

A different but related issue is connected to the obser- 
vation that indeed, the set of entangled states is dense 
in the set of all states for continuous variable systems 
[23^ Hence one may ask how realistic the separable 
states found here are or whether they are artificial math- 
ematical constructs. In other words, the actual state of 
the system may never be the exact thermal state, which 
would be separable above T C rit, but an entangled non- 
equilibrium state nearby. This question is identical to 
asking how clear a cut exists between a quantum and 
classical world. The paper [23[ has a very interesting 
point to make about this: namely, that as long as we 
only look at states of bounded energy (say, of the same 
scale as the thermal state under consideration) in the 
vicinity of the thermal state, then the amount of entan- 
glement in theses states is bounded by a function of the 
distance (measured via the trace norm). In other words, 
while we cannot be sure that the actual state, fluctuating 
a tiny bit away from the exact thermal state, is indeed 
separable, we do know that the amount of entanglement 
in it must be very small - as long as we can control its 
mean energy and distance from the thermal state. 

Finally, the choice of the 'parties' to find entanglement 
between, which we addressed in Remark I, becomes even 
more urgent for the continuum limit of the harmonic 
model. In the limit of infinitely many harmonic oscil- 
lators each site turns into a point in space that can os- 
cillate and the system becomes a quantum field. The 
entanglement can now exist between modes of the field 
instead of single particles [2(| [27| and one is left with the 
choice of which modes to talk about (2^|. The decou- 
pled, unentangled normal modes in the lattice, the lat- 
tice vibrations, become continuous single-particle encrgy- 
cigenfunctions and the entanglement between the former 
discrete sites in the harmonic chain becomes entangle- 
ment between points in space in the quantum field. In 
a rescaled and non-diverging theory the interaction be- 
tween spatial points, w, is infinitcsimally small and hence 
no entanglement can be found between them. However, 
analogous to the block entanglement in finite harmonic 
chains, discussed in [20| . the entanglement between finite 
volumes in space can be tested as done in 0, HI1 • 
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